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Abstract 

Opposing to a (common) belief against the existence of a thermodynamic-like potential for the 
KPZ equation, here we present a derivation for such a functional. With its knowledge we prove 
some global shift invariance properties previously conjectured by other authors. The procedure 
could be extended in order to derive a more general form of such a functional leading to other 
known related nonlinear kinetic equations. Exploiting the KPZ's functional, and for arbitrary 
dimension, we have obtained the exact form of the stationary probability distribution function 
and have shown a couple of examples of how it is possible to exploit it in order to obtain relevant 
results like finding support to the conjecture that in the strong coupling regime a critical dimension 
doesn't exists. 
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Phenomena far from equilibrium are ubiquitous in nature, including among many other, 
turbulence in fluids, interface and growth problems, chemical reactions, biological systems, 
as well as economical and sociological structures. During the last decades the focus on 
statistical physics research has shifted towards the study of such systems. Among those 
studies, the understanding of growing kinetics at a microscopic as well as on a mesoscopic 
level constitutes a major challenge in physics and material science . Some recent 

papers have shown how the methods and know-how from static critical phenomena have 
been exploited within nonequilibrium phenomena of growing interfaces, obtaining scaling 
properties, symmetries, morphology of pattern formation in a driven state, etc 0,00, Sf 



i4 

It is a common belief that the nontrivial spatial-temporal behavior occurring in several 



nonequilibrium systems 



111 ], originates from the non-potential (or non-variational) char- 



acter of the dynamics, meaning that there is no Lyapunov functional for the dynamics. 
However, Graham and co-workers have shown in a series of papers [3] that a Lyapunov-like 
functional exists for very general dynamical systems, like the complex Ginzburg-Landau 
equation. Such a functional is formally defined as the solution of a Hamilton- Jacob i-like 



equation, or obtained in a small gradient expansion |12l . [13]. The confusion associated 
with the qualification of nonvariational dynamics comes from the idea that the dynamics 
of systems having nontrivial attractors (limit cycle, chaotic) cannot be deduced from the 
minimization of a potential playing the same role as the free energy in equilibrium systems 



16]. Nevertheless, this does not preclude the existence of a Lyapunov functional for 



the dynamics that will have local minima identifying the attractors of the system. How- 
ever, once the system has reached an attractor that is not a fixed point, the dynamics 
proceeds inside the attractor driven by nonvariational contributions to the dynamical flow, 
that do not change the value of the Lyapunov functional, implying that the dynamics is 
not completely determined once the indicated functional is known. This situation has some 
known examples even in equilibrium statistical mechanics 13]. Hence, the Lyapunov func- 
tional, or nonequilibrium potential (NEP) 121 ]. plays the role in nonequilibrium situations 
of a thermodynamic-like potential characterizing the global properties of the dynamics: at- 
tractors, relative (or nonlinear) stability of the attractors, height of the barriers separating 
attractions basins, and offers the possibility of studying transitions among the attractors 
due to the effect of (thermal) fluctuations. 
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In a recent series of papers we have shown several results related to the obtention and ex- 
ploitation of the indicated NEP's concept in scalar and non-scalar reaction-diffusion systems 
(see [l6] and references therein). In particular we have exp loited those results for the study 
of stochastic resonance [l9] in extended systems (see 15, 16, 20] and references therein). 
In those works, we have analyzed problems of stochastic resonance in scalar and activator- 
inhibitor systems, systems with local and nonlocal interactions, system-size stochastic reso- 
nance, etc. 

Here, and related to the kinetics of growing interfaces, we discuss the case of the Kardar- 
Parisi-Zhang equation (KPZ) 2jJ. This equation describes the evolution of a field h(x,t), 
that corresponds to the height of a fluctuating interface, 



d A 

— h(x, t) = zA7 2 /i(x, t) + - (Vh{x, t)f + K Q + £(x, t), 



(1) 



where £(x,t) is a Gaussian white noise, of zero mean t)) = 0) and correlation 

t)£(x', t')) = 2e5(x — x')5(t — £'). As indicated above, this nonlinear differential equa- 
tion describes fluctuations of a growing interface with a surface tension given by u, A is 
proportional to the average growth velocity and arises because the surface slope is parallel 
transported in such a growth process. 

Opposing to a claim in a recent paper The KPZ equation is in fact a genuine kinetic 
equation describing a nonequilibrium process in the sense that the drift vV 2 h + f (V/i) 2 ~ F 
cannot be derived from an effective free energy] we show here that such a nonequilibrium 
thermodynamic-like functional (NETLP) for the KPZ equation exists. Exploiting its knowl- 
edge, we will discuss conjectures advanced in [5| and how they are fulfilled. We also briefly 
discuss how to extent the derivation procedure in order to consider more general forms of 
kinetic equations. Finally, we obtain the stationary (or asymptotic) probability distribu- 
tion function (pdf) -valid for any dimension and unknown till now- and also derive the 
form of the NEP and, through the analysis of a couple of simple examples, we discuss how 
the knowledge of this NETLP and pdf could be exploited in order to obtain some relevant 
results. 

The Lyapunov functional or NETLP for the KPZ equation is given by 



4z/ 



-K + ^(Vh(x,t)f 



dx. 



(2) 
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It is easy to prove that this functional fulfills both, the relation 

as well as the (Lyapunov) property ^^[h] = —T[h] {j^^h]) 2 < 0, where the function T[h] 
is given by 

^Ve-tM^). (4) 

Hence, from this free energy-like functional, and by a functional derivative, we can obtain 
the KPZ kinetic equation. It corresponds to a relaxation model, analogous to model A 
according to the classification in Hohenberg & Halperin's review [171 ] . 

It is worth to make here a remark. In one hand, in the standard "model A" it is known 
that the dynamics can be seen as a superposition of modes that decay exponentially towards 
a steady state, with time-dependent correlations obeying some constraints such as positivity. 
But on the other hand, in the KPZ problem it is known that the relaxation of perturbations 



decay following a stretched exponential form 



22 
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In order to show how to obtain the above indicated functional, we start considering the 
following simple scalar reaction-diffusion equation for a positive (<fi > 0) field (p(x,t), as it 
corresponds to a probability density, 

d 

—<f){x, t) = z/V 2 0(x, t) + a 4>(x, t) + r)(x, t) (f)(x, t), (5) 

where a is a constant, and rj(x, t) is also a Gaussian white noise of zero mean, and intensity 
a, and we assume the Stratonovich interpretation. It is well known that the system in Eq. 



has the following NETLP 



26] 



*oM= / {-^(x,t) 2 + ^(V^x,t)) 2 ]dx, (6) 



where Q indicates the integration range. As has been shown in previous works 

d_ 

at- 



in 



addition to fulfilling the Lyapunov propertiy ^■^ r [0] < 0, it also fulfills the relation 



where the contribution from the boundaries is null, due to the variation 5(f> being fixed there 
(= 0), as usual. 
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Let us now define a new field, h(x, t), that as indicated before corresponds to the interface 
height, exploiting the so called Hopf-Cole transformation h(x,t) = ^ln0(x,t), with the 
inverse <f)(x,t) = e^ h ^ x,t \ It is straightforward to show that exploiting this transformation, 
the original Eq. (jSJ) becomes Eq. (CQ), with a = K j- and a = ^e. However, the noise term 
that in Eq. (jSJ) has a multiplicative character, in the transformed Eq. (CQ) becomes additive. 

If we now apply the same transformation to the NETLP indicated in Eq.©, it is immedi- 
ate to obtain the functional shown in Eq. (j2J). Hence we have a free energy-like functional 
from where the KPZ kinetic equation can be obtained through functional derivation. Clearly, 
the contribution to the variation that come from the boundaries is again null. 

It is worth to consider once more Eq. (J5j), but now including a typical limiting term of 
the form: —b(p(x,t) 3 . The resulting reaction-diffusion equation corresponds to a version of 
the so called Schlogl model The associated NETLP will have an extra term of the form 
+| J n (j)(x,t) 4 dx. Applying once more the previously indicated Hopf-Cole transformation, 
in Eq. (CD) a new associated term arises, having the form —7 e«^f^ (b = 7^;). The new 



equation corresponds to a form of the so called bounded-KPZ [28|, [29J]. Clearly, we will 
also have an extra term in the associated NETLP (Eq. (j2J)). However, in what follows we 
consider the case b = 0, analyzing only the more "usual" form of the KPZ equation indicated 
byEq. ©. 

Let us now check some of the properties previously assumed for such a functional. Ac- 
cording to the analysis of global shift invariance in 5J, it is easy to see that the relations 
indicated by Eq. (9) in p] are fulfilled. That is, we can readily prove that if I is an arbitrary 
(constant) shift 

F[h + l]=K[l]F[h]; T[h + l} = K[l}- 1 V[h], (8) 

with K[l] = e^(~ Til]- 1 ). 

To prove other conjectures also indicated in [5], we introduce the free energy-like density 
^[hj'Vh], which is defined by T[K\ = J dxJ-[h,Vh]. The relations we refer are 

T[h,Vh] = e^![(V/i) 2 ]; T[h,Vh] = e^T^Vh) 2 ]. (9) 

According to the form of the NETLP indicated in Eq. and the definition of J-[h, Vft], 
it is clear that the first relation above results "trivially" true. For the second relation we 
have that T[h, V/i] = e~ sh<yX ^T , where r o = 1, and s = -, as T[h] is not function of V/i. In 



addition, it can be also proved that the indicated NETLP is invariant under the nonlinear 

Q 

Galilei transformation that, as discussed in [7], are fulfilled by the KPZ equation. 

We can go still further and look for the possibility of deriving the NETLP for other 
forms of related kinetic equations. To do that, we should assume that we have a non-local 
react ion- diffusion equation, as in 

d 



dt 



<p(x,t) = uV 2 (p(x,t) + a(f)(x,t) - j3 / dx'G(x,x')(j)(x',t) + cf)(x',t)ri(x',t), (10) 



where, as discussed in 20[ the kernel G(x, x') could be of a very general character, and (3 is 
the interaction intensity. 

As we have done before, using the Hopf-Cole transformation we could obtain a generalized 
{nonlocal) form of the KPZ equation. However, in order to obtain sensible results we should 
assume that the translational invariant kernel (G(x, x') = G(x — x')) is of a very short range 
and has some symmetry properties. What results is a nonlocal contribution that, even 
differing from the one discussed in 3jJ, l32|, is of much interest and relevance. Repeating 



the previous procedure we find the associated NETLP, but we will not discuss this aspect 
here (see js(3]). We only want to remark that the contributions that arise have the same 
form of those ones that arose in several previous works, where scaling properties, symmetry 
arguments, etc, have been used to discuss the possible contributions to a general form of 



the kinetic equation 



33[. The different contributions that arise are tightly related 



34|, the 



to several of other previously studied equations, like the Kuramoto-Sivashinsky 
Sun-Guo-Grant (SGG) equation [3^], and others as in |s| and 10]. 

The knowledge of the NETLP for the KPZ equation allows us to readily write the asymp- 
totic long time probability distribution function (pdf), valid for any dimension, which (due 
to the "diagonal" character of T[h]) is given by 



V as [h(x,t)) ~ exp ^ - 



exp 



exp < - 



h(x,t) 



dx 



dx 



h r 



A 



dijj uV 2 i/j(x, t) + - (W>(5, t)) 



h r Sijj 



T[h\F[h] 

e, V \-^-jdx{Vhf + ±jdxj h 



h(x,t) 



d^j 



exp 



$[h) 



(11) 



where, in order to simplify, we assume K = 0, and h r is an arbitrary reference state. The 
second line results from just replacing the relation in Eq. (JS]), while the third and fourth lines 
results transforming the pdf for the original field, or by using functional methods. The first 
part of the last line shows a nice structure, where we can identify a contribution, Gaussian on 
the slope, plus a "correction" proportional to A. For the one-dimensional case, it is possible 
to show that only the well known Gaussian result survives [2|, 13J , if the adequate boundary 
conditions are taken into account. Otherwise, the indicated expression is the complete 
solution for the pdf, irrespective of the boundary conditions. 
The last line gives us the definition for that is 

^ f fh(x,t) 



$[h] = ^jdx (Vh) 2 - '- J dx J dip {Vipf , 

a functional that could be identified as the nonequilibrium potential (NEP) lil 
fulfills 

6$[h] 



(12) 



13 



15| that 



9 w v 



8h(x, t) 



(13) 



as well as -§i&[h] 



tion of other techniques 12 



m 



h] ) . The knowledge of $[/t] opens the door to the exploita- 



151 ] . offering alternatives to the usual renormalization group 



ones. Clearly, it also shows that the claim done in the phrase indicated at the beginning [3] 
is not correct. 

In order to show the possibilities that offers the knowledge of such a NETLP, we now 
present a simple example. Let us consider a sightly different case, where we only have a 
spatial quenched noise ( "disorder" ) instead of the spatio-temporal noise we have considered 
so far. The associated equation is 



d X 

—h{x, t) = vV 2 h(x, t) + - {Vh(x, t)) 2 + K Q + 



(14) 



where, as in previous studies 
For this case we have that 



36 







371 . |38| . is a quenched, Gaussian distributed, noise. 



h(x, t) 



-v[h\- 



dt ' 5h(x, t) ' 

with the same form of as in Eq. (T5]), but K a replaced by K Q + -d(x) 



(15) 



The indicated studies have shown that the front profile presents a triangular structure 



38j | and, clearly, it is of relevance to determine the slope of such structures. Using the known 
form of the NETLP, we can minimize this free-energy-like functional and obtain, in the 1-d 

1 /3 

case, that such a slope is a = (^4) , a value that agrees quite well with the numerical 
evaluations. 

Summarizing, we have here found the form of the Lyapunov functional or NETLP for the 
KPZ equation. From this NETLP, and through a functional derivative, we have obtained 
the KPZ kinetic equation, and have also shown that it fulfills global shift properties, as 
well as other ones anticipated for such an unknown functional. Even more, it is possible to 
extend the procedure to derive such a functional, considering nonlocal terms, and in such a 
way derive more general forms, that includes several other kinetic equations studied in the 
literature of interface growing phenomena. 

As indicated in the literature, dynamic renormalization group techniques, being useful 
and powerful, in many cases only offers incomplete results, having no access to the strong 



coupling phase |2|, [39] . Hence, it is clear the need of alternative, complementary, ways 



to analyze the KPZ and related problems (as an example see [32, |40( for a self-consistent 
expansion). The present results clearly open new possibilities of making non-perturbational 
studies for the KPZ problem. We expect that through the analysis of long time mean values 
of h(x,t), or of correlations, we could extract information about scaling exponents. Such 
study will be the subject of forthcoming work |4l| ] . 
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